Wormhole solutions in the Randall- Sundrum scenario 

M. La Camera]] 

Department of Physics and INFN - University of Genoa 
Via Dodecaneso 33, 16146 Genova, Italy 

Abstract 

In the simplest form of the Randall-Sundrum model, we consider the metric generated by a 
static, spherically symmetric distribution of matter on the physical brane. The solution to the 
five-dimensional Einstein equations, obtained numerically, describes a wormhole geometry. 

Keywords: Brane world, wormholes. 



I. INTRODUCTION 



In the last few years large extra dimensions and brane worlds have been the subject 
of intensive investigations. The situation can be simplified to a five-dimensional problem 
where matter fields are confined to the four- dimensional spacetime while gravity acts in 
five dimensions. If the extra dimension has a finite extension its compactification radius 
may be larger than the Planck length without conflict with observations [1]. Alternatively, 
the extra dimension may be kept warped and compactified or uncompactified [2,3] or even 
neither compact nor warped with a truly infinite size [4]. A considerable effort has been 
devoted in exploring possible phenomenological and observable consequences of the brane 
world scenario. In particular Lorentzian wormholes [5,6], smooth bridges connecting different 
universes or remote parts of the same universe, have gained much attention and have found 
a natural source in brane world models (see [7] and references quoted therein). In this letter 
we will show that a solution to Einstein equations in the Randall- Sundrum scenario with a 
mass point on the physical brane can describe a wormhole geometry. The paper is organized 
as follows. In Section II we start from the simplest form of the Randall-Sundrum model 
and modify the corresponding metric, which is asymptotically anti-de Sitter, to a metric 
generated by a static, spherically symmetric distribution of matter on the physical brane or, 
equivalently, axially symmetric in the bulk. By a suitable choice of coordinates we succeed 
in reducing the solution of the five-dimensional Einstein equations to the solution of a single 
ordinary differential equation in one unknown function. In Section III we solve numerically 
the highly nonlinear equation obtained and show that there are wormhole solutions in the 
brane world. Some concluding remarks are given in Section IV. 

II. THE MODEL 

There are in fact two models due to Randall and Sundrum with essentially the same 
framework. In the RSI model [2] we have two 3-branes with equal and opposite tensions 
which rigidly reside at the boundaries (y — and y = r c n) of a slab of a five- dimensional 
anti-de Sitter bulk spacetime of radius £. The five- dimensional Einstein equations in the 
vacuum are 

Gab = Rab - ~9abR = - ~^jz^ AB W 
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Here M 5 and A denote respectively the five- dimensional Planck scale and the (negative) bulk 
cosmo logical constant. The Randall- Sundrum solution to the previous equations is 

ds 2 = e~ 2k lyl ri fMU dx»dx u + dy 2 (2) 

where r)^ is the flat Minkowskian metric, and k — \ is a scale factor for the extra dimension 
y. The above solution holds only if the tensions of the branes and the bulk cosmological 
term are fine-tuned and in particular this implies that k must take the value k = 
Moreover the following relation is derived 

M 2 A = ^L(\ - e ~ 2kr ^) (3) 

which gives a well defined value for the 4D Planck scale M4 even in the infinite radius limit. 
So if M 5 and k are taken both of order M 4 and kr c ii xs 35 the model suggests a solution to 
the hierarchy problem by showimg that TeV energy scales on the "TeV brane" at y — r c ir 
correspond to 10 16 TeV energy scales on the "Planck brane" at y — 0. In an alternate 
scenario known as RS2 [3] the second brane, where now is supposed to be concentrated 
the graviton zero mode, is taken off at infinity letting r c — > 00, while the Standard Model 
particles are assumed to live on the brane at y — 0. Such localization of gravity can obviate 
the need for the compactification of the fifth dimension. RS2 does not solve the hierarchy 
problem as RSI does, but it is of interest from a purely gravitational point of view. While 
in the absence of matter the line interval can be exactly obtained, complications arise in the 
presence of matter fields, and even if an analytic expression could be found for the metric 
produced by a matter distribution on the physical brane it certainly would involve special 
functions, as one can see in linear and beyond linear order approximate solutions [8,9,10]. 
Here we assume for simplicity that the matter fields reduce to a point mass on the physical 
brane. Having in mind the four-dimensional McVittie solution [11], which represents the 
gravitational field of a mass particle in a Friedmann-Robertson- Walker universe and its 
generalization to higher dimensions [12], we employ similar methods to obtain a solution 
although the metric we are searching for is not spherically symmetric in five dimensions. 
We start from the sourceless AdS$ solution given by 

ds 2 = e~ 2k ^ {dr 2 + rW - dt 2 ) + dy 2 (4) 

where dfl 2 is the line element on a two-dimensional sphere of unit radius, and apply to it 
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the following transformation of the Robertson type [13] 

P 



kz 



y/l + k 2 p 2 



2 - — log (1 + k 2 p 2 ) 



(5) 



so obtaining 



ds 2 = 



dp 2 



p 2 dn 2 - e~ 2kz (1 + k 2 p 2 ) dt 2 + (1 + k 2 p 2 ) dz 



(6) 



1 + k 2 p 2 

In the presence of a point mass m bound at r = to the physical brane in the old coordinates 
and consequently bound at p = in the new coordinates, we make the ansatz that the line 
element (6) becomes 



ds 2 = e a ^ dp 2 + p 2 d£l 2 - e 2kz+ ^ dt 2 + e<"W dz 2 



(7) 



In the above five- dimensional metric we used only four of the five available degrees of freedom, 
so (7) still has a degree of generality which can applied in future contexts. From Einstein 
equations (1), written in the mixed form 



G B A = Qk 2 5 B A 



we obtain, omitting for simplicity the p-dependence 

G? = !e w (V-u/) = 



(8) 



(9) 



Therefore putting 7 = 0; one realizes by inspection that two of the components of the mixed 
Einstein tensor, namely G\ and G z z , are equal. 
Finally the relevant equations are 
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Ap 



2 { e - a [A - 4e a ~P + 4e Q - 7 AV + p 2 (3' 2 + 8r 7 ' + p 2 7 ' 2 
+4p/?(l + P7 / )]} = Qk 2 



—{e- a [Ae a ^k 2 p + pf3' 2 + 4 7 ' + 3p 7 ' 2 + 2/7(2 + py') 
4p 

-a' (2 + pp + 2p 7 ') + 2pP" + 4p 7 "]} = 6k 2 

-^r{e- a [-4e a -P + 4 + 3p 2 f3' 2 + 4p 7 ' + p 2 7 ' 2 + 2p/3'(6 + pi) 
Ap z 

-pa'(A + 2pp + pi) + Ap 2 {3" + 2p 2 7 "]} = Qk 2 



(10a) 
(10b) 
(10c) 



Solving for e a in equation (10a) and substituting this value into equations (10b) and (10c) 
we obtain: 

a _ ^ [(2 + P_P) (2 + p? + Wj + pV 2 ] 
4 [e^ + e/ 3 (-l + 6eT) Fp 2 ] 

3V 72T$(P) = (lib) 



where 



[(2 + pp) (2 + + 4p 7 ') + p 2 y 2 ] 

3 ( 2 + ^) m = (He) 

[(2 + p (2 + p p> + 4p y ) + p 2 y 2 ] yi > y J 

$(p) = e^p{i[p[p(3' 2 + 4(p 7 ' + 1)P' + i {pi + 8) - 2pP"] + 8] 
+2p(pP> + 2) 7 "} - Pp 2 {4[p[p 2 /?' 3 + 3p(p 7 ' + 2)f 
-miiPl' - 4) + P7'1 - 4]/?' - iWiPl' + 6) - 3p/3"] - 6] 
-6p 7 "] + 8] - e^[p[p 2 /5 /3 + 2p(2p 7 ' + 3)/3 /2 + [p^^pV + 16) 
-2p 7 "] + 12]P' + 2 7 '[p( 7 ' + pP") + 6] - 4p 7 "] + 8]} (12) 

Equation (lib) can be satisfied by 7 = 0, but in this case (11c) and (11a) become respectively 
k 2 = and e a = e^{2 + pf3') 2 so it is immediate to verify, passing to standard coordinates, 
that we would have a flat five-dimensional space solution without cosmo logical constant. In 
the same way equation (11c) can be satisfied by P = 21n(p /p) but, apart the unpleasant 
fact that the standard radial coordinate should have the constant value p , (lib) becomes 
1 + po 2 ^ 2 = so cannot be satisfied. Finally the choice 7 = and P = 21n(p /p) is to be 
ruled out because it makes not soluble already the system (10) from which we started. As 
a consequence (lib) and (11c) are both satisfied only by $(p) = 0. The fact that now we 
are faced with the system (11a) and (12) of two equations in the three unknown functions 
a,P,j is not surprising because in writing the line element (7) we retained the possibility, 
due to the coordinate freedom of Einstein equations, of still applying one degree of freedom 
so we can arbitrarily reduce the unknown quantities from three to two. A possible choice 
can be made exploiting the solutions to the system in some limiting cases. If m = and 
k > the solution is given by the line element (6), while if m > and k = the solution has 
been investigated by various authors [14-17] and the corresponding line element, written in 
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isotropic form, is given by 



ds 2 = 1 + 



\/3' 



m 



2p 



V3 



m 



2p 



1 + 



rn 



2p 



2(l-2/ v / 3) 



(dp 2 + pW) 



1 - 



v 7 ^ 



rn 



2p 



1 + 



>/3: 



rn 



2p 



2/V3 



dt 2 



+ 



V orr 



1 + 



v/3 



m 



2p 



*/ 2 



(13) 



The above line element was obtained in [17] as a solution to the five- dimensional Einstein 
equations outside a static and three-dimensional spherically symmetric distribution of mass 
m and scalar charge a in equal amount and there it has been shown that it describes a 
Lorentzian wormhole with throat at p = \/3/2m. Now comparing the p-dependence of the 
coefficient multiplying dt 2 in the line elements (6) and (13), where respectively m = and 
k — 0, a suitable but of course not the only one choice for the function e 7 ^ when both m 
and k are different from zero may be the following: 

2/V3 



1 - 



rn 



2p 



1 + 



rn 



2p 



;i + k 2 P 2 ) 



(14) 



So equation $(p) = reduces to a second order differential equation, which for simplicity 
is not written here, in the unknown function f3(p) and will be solved numerically. As to 
the boundary conditions for j3 and its first derivative, they will be fixed once we know the 
behavior of (5(p) when mj p is sufficiently smaller than unity and both m and k are different 
from zero. This can be obtained from the equation $(p) = inserting in it the expansion 
of 7(p) in those conditions and solving by series the resulting differential equation for f3(p). 
We find, with the requirement that (5 vanishes at infinity 



2m \ 2 fc 2 J 

P P A 



(15) 



Opposite to the case (m > 0, k — 0), where (5 = 4m/ p + (3/2m 2 )/p 2 + • • • , now (5 reaches 
its asymptotic value (3 = from below. 
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III. NUMERICAL CALCULATIONS 



In view of searching for a numerical solution we have to fix numerical values for the two 
parameters k and m. Now in both the models of Randall-Sundrum one can derive from 
equation (3) the following relation between k and the fundamental mass scales 

Mf v ' 

so with M 5 and k taken, as in RSI, of order M 4 10 19 GeV, one obtains k ps 10 32 mm' 1 

and since kr c n p=* 35 it follows that r c p=* 10~ 31 mm, so this is not a model with a large 

extra dimension. There are however various extensions of the generalized RS2 model (see 

e.g. [18] and references quoted therein) where the extra dimension is compact and much 

larger than Planck scale. In any case experiments impose the bound r c < 1 mm, which 

implies k > 10 mm -1 , so we shall let k to vary by many orders of magnitude. In searching 

for the numerical solution to equation <3>(p) = we will pass to the adimensional variable 

p/m and as a consequence in that equation m and k will appear only in the adimensional 

quantity mk, but clearly not all the possible combinations of m and k leading to the same 

value of mk are physically acceptable. As a numerical example if m refers to the Sun mass 

then m ps 10 6 mm and therefore 10 7 < mk < 10 38 . The precision of calculations were tested 

verifying the reproduced value of the quantity 6k 2 in the right-hand side members of Einstein 

equations (10b) and (10c) in a wide range of the variable p. The relative error proved to 

be a decreasing function of p, some of its approximate values being 10~ 2 at p ^ m, 10~ 4 at 

p = 2m and 10~ 6 at p = 10 6 m. 

To discuss the numerical solution of equation (12) let us consider the plots in figure 1. In 
the case m > and k = the function e@, given now by equation (13), goes to infinity as p 
goes to v / 3/2m and decreases to one as p — > oo. If m = and k > one has e 13 = 1 from 
equation (6). Finally when both m and k are different from zero e 13 becomes again infinite 
at p = \/3/2m, reaches a minimum less than unity near p = 1.15 m then increases towards 
the asymptotical value e 13 = 1. It is worth noticing that the last behavior is numerically 
the same, in the limit of machine-precision, for all the values of mk in the above considered 
range and this happens because the values mk ^> 1 make equation $(p) =0 very weakly 
dependent on mk. As to the function e a given by equation (11a) we found, using the 
calculated values of (5 and for all values of mk as before, that it is positive when p is greater 
than p min p^ (y/3/2+e) m, where e ps 2.5 10 -4 . Our solution is therefore defined for p > p m in- 



m=0, k>0 



m=0, k>0 
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FIG. 1: The function for different FIG. 2: The function S for different 
values of the pair (m,k). values of the pair (m,k). 

At Pmin the functions e 13 and e 7 are finite, different from zero and both much greater than 
unity. Finally, for p min < p < oo the function e a results decreasing and of order l/(k 2 p 2 ), 
except in the interval (0.95 m, 1.12 m) where it is slightly increasing. To see what is going on, 
look at the standard radial coordinate R(p) = pe 13 ^ 2 depicted in figure 2 for different values 
of the pair (m, k). The case mk > lwe are treating here is similar to the case (m > 0, k = 0) 
of equation (13) already discussed in [17] and will bring to analogous conclusions. More in 
detail, consider the proper circumference of radius p given by C(p) = 2tcR(p). We have 
that C(p) becomes infinitely large as p — > oo and asymptotically very large, mantaining 
however finite, as p — > p m in and has a minimum at po ~ 1.04 m. This numerical value was 
estimated plotting the derivative dR/dp versus p and looking at the value of p where it 
vanishes. The regions p G (p m in, Po) and p G (po,oo) are two spacetimes associated with a 
wormhole whose throat occurs at p where C(p) is minimum. The wormhole is in this case 
aymmetric under the interchange of the asymptotic surfaces (p = p m in and p = oo) and 
traversable in principle due to the absence of event horizons. We verified numerically that 
the Kretchmann invariant becomes infinity only at p = i/3/2m which is not part of the 
manifold. 
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IV. CONCLUSION 



Some remarks seem here to be appropriate. First, it is known that a wormhole geometry 
can appear only if some of the energy conditions are violated and this requires the presence 
of some amount of "exotic" matter. The violation could be justified incorporating the cos- 
mological constant into an energy-momentum tensor describing the vacuum contribution to 
the density and pressure. However wormhole solutions can be obtained also with vanishing 
cosmological constant in a Ricci flat spacetime, where energy conditions are trivially satis- 
fied, of dimension iV > 5 as in [17] where it was shown that every time an extra dimension is 
reduced a massless scalar field appears and plays the role of "exotic" matter. Also we would 
like to stress that in the framework of induced matter theory [19] the four-dimensional Ein- 
stein equations with sources can be locally embedded in five-dimensional Einstein equations 
without sources so although the matter seems exotic in four dimensions the five-dimensional 
spacetime on both sides of the throat is that of a vacuum and energy conditions are again 
trivially satisfied. Second, according to the common assumption made both in brane world 
and in induced matter theories, our five-dimensional metric tensor depends explicitly on the 
fifth coordinate, as can best seen inverting transformation (5) and applying it to the line 
element (7). One obtains 

ds 2 = e- 2fc W { -JJLIJ—— }dr 2 + e^dil 2 - -^———dt 2 } 

U (1 + k 2 r 2 e~ 2k \y\) 2 \ (1 + k 2 r 2 e- 2k \y\) J 

eWe-*' + (1 + ,,X_ 2%I)2 ] V + 2kre~^ [ - ? 4- (l - pgpg^ ] driy (17) 



+ 



where a, (5 and 7 denote the functions a, (3 and 7 written in the coordinates r and y. 
In this short letter we have not considered more general gauge transformations. Further 
investigations should have to take into account, as pointed out in [20], the influence that 
coordinate transformations may have on the interpretation of physical quantities in one 
lower dimension when in the reduced dimension is still present the extra coordinate. 
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